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1.  Introduction 

An  n- variable  Boolean  function  /  is  a  map  from  the  n  dimensional  vector 
space  =  {0,  l}n  into  the  two-element  field  F2,  that  is,  an  n-variable 
Boolean  function  /  is  a  multivariate  polynomial  over  F2.  Denoting  the 
addition  operator  over  F2  by  '+’,  a  Boolean  function  can  be  thought  as  a 
multivariate  polynomial,  called  the  algebraic  normal  form  (ANF) 

f  (x  1,  •  •  •  ,  X,f)  —  ao  4“  ^  '  CtiXi  ^  '  CtijXiXj  +  •  •  •  +  Q-12.. ,nX\X2  •  •  •  Xn, 
l<i<n  1  <i<j<n 

where  the  coefficients  ao,  ay, . . . ,  ai2...n  £  F2.  The  maximum  number  of 
variables  in  a  monomial  is  called  the  ( algebraic )  degree,  and  it  is  denoted 
by  deg(f).  If  all  monomials  in  its  ANF  have  the  same  degree,  the  Boolean 
function  is  said  to  be  homogeneous. 

Functions  of  degree  at  most  one  are  called  affine  functions.  An  affine 
function  with  constant  term  equal  to  zero  is  called  a  linear  function.  Define 
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the  scalar  product  of  x  =  (aq,  , . . ,  xn),  y  =  (y\, . . .  ,yn)  both  in  ¥%,  by 
x  •  y  =  E?=1  xlyl.  The  (Hamming)  weight,  denoted  by  wt(x.),  of  a  binary 
string  x  is  the  number  of  ones  in  x,  and  the  Hamming  distance  c?(x,  y) 
between  x  and  y  is  the  number  of  positions  where  x,  y  differ.  An  n- variable 
function  /  is  said  to  be  balanced  if  its  output  column  in  the  truth  table 
contains  equal  number  of  0’s  and  l’s  (i.e.,  wt(f )  =  2n_1).  The  nonlinearity 
of  an  n- variable  function  /  is  the  minimum  distance  to  the  entire  set  of  affine 
functions,  which  is  known  to  be  bounded  from  above  by  2re-1  —  2n/2~1. 

We  define  the  (right)  rotation  operator  pn  on  a  vector  (x\,x2,  ■  ■  ■ ,  xn)  € 
F2  by  Pn(x i,X2,  ■  ■  ■  ,xn)  =  (xn,xi,x2,  ■  ■ .  ,xn-i).  Hence,  p\  acts  as  a  k- 
cyclic  rotation  on  an  n-bit  vector.  A  Boolean  function  /  is  called  rota¬ 
tion  symmetric  if  for  each  input  (x\, . . .  ,xn)  in  ¥%,  f(plf(x  ±,...,xn))  = 
f(x  1, . . . ,  xn),  for  1  <  k  <  n.  That  is,  the  rotation  symmetric  Boolean  func¬ 
tions  (RSBF)  are  invariant  under  cyclic  rotation  of  inputs.  A  partition  of 
some  cardinality  gn  is  generated  by  Gn(x  1, . . .  ,xn)  =  {p^{x  1,  •  •  •  ,xn)\l  < 
k  <  n},  and  so,  the  number  of  n- variable  RSBFs  is  29n .  It  was  shown  [11] 
that  gn  =  ^  Efc|n  </*(&)  2* ,  where  <f  is  Euler’s  totient  function.  We  refer 
to  [8,  9,  11]  for  the  formula  on  how  to  calculate  the  number  of  partitions 
with  weight  w,  say  gn>w,  for  arbitrary  n  and  w. 

A  rotation  symmetric  function  f(x  1, . . . ,  xn)  can  be  written  as 

Uo  0  1  V-"  1  T  ^  ]  Cl\jX\Xj  -I-  -t- 

where  ao,  a\,  a±j, . . . ,  ai2...n  £  ^2j  and  the  existence  of  a  representative  term 
x\ Xi2  . . .  xH  implies  the  existence  of  all  the  terms  from  Gn(x\Xi2  . . .  Xit )  in 
the  ANF.  This  representation  of  /  (not  unique,  since  one  can  choose  any 
representative  in  Gn(xixl2 . . .  x^))  is  called  the  short  algebraic  normal  form 
(SANF)  of  /.  If  the  SANF  of  /  contains  only  one  term,  we  call  such  a 
function  a  monomial  rotation  symmetric  (MRS)  function.  Certainly,  the 
number  of  terms  in  the  ANF  of  a  monomial  rotation  symmetric  function  is 
a  divisor  of  n  (see  [11]). 

We  say  that  two  Boolean  functions  /(x)  and  g(x)  in  Bn  are  affine  equiv¬ 
alent  if  g(x)  =  /(xA  +  b),  where  A  G  GLn(¥2)  (n  x  n  nonsingular  matrices 
over  the  finite  field  F2  with  the  usual  operations)  and  b  is  an  n- vector  over 
F2.  We  say  /(xA  +  b)  is  a  nonsingular  affine  transformation  of  /(x).  It  is 
easy  to  see  that  if  /  and  g  are  affine  equivalent,  then  they  have  the  same 
weight  and  nonlinearity:  wt(f)  =  wt(g )  and  Nf  =  Ng  (these  are  examples 
of  affine  invariants) . 

There  are  cases,  when  it  is  known  that  these  invariants  are  also  sufficient 
(two  quadratic  functions  are  affine  equivalent  if  and  only  if  their  weights 
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and  nonlinearity  are  the  same-see  [3],  for  example).  However,  in  general, 
for  higher  degrees,  that  it  is  not  the  case. 


2.  Background  on  S'-equivalence 

In  [2]  the  authors  introduced  the  notion  of  S -equivalence  f  ~  g,  which 
is  the  affine  equivalence  of  monomial  rotation  symmetric  (MRS)  functions 
/,  g  under  permutation  of  variables  (we  will  write  here  /  ~  g,  for  easy 
displaying) . 

An  n  x  n  matrix  C  is  circulant ,  denoted  by  C(ci,C2, . . .  ,cn),  if  all  its 
rows  are  successive  circular  rotations  of  the  first  row,  that  is, 


( C1 

C2  • 

•  •  cn  \ 

c  = 

Cn 

Cl  . 

•  •  cn— 1 

\C2 

C3  • 

..  d  ) 

On  the  set  Cn  of  circulant  matrices  an  equivalence  relation  was  introduced 
in  [2]:  for  A\  =  C(ai, . . . ,  an),  =  C(&i, . . . ,  bn),  then  A\  ps  A^  if  and  only 
if  (ai, . . . ,  an )  =  Pn(bi, . . . ,  bn),  for  some  0  <  k  <  n  —  1.  It  was  shown 
that  the  set  of  equivalence  classes  (the  equivalence  class  of  C(ai,  02, . . . ,  an) 
is  denoted  by  C(ai,  a%, . . . ,  an),  or  (C(ai,  02, ... ,  an)))  form  a  commutative 
monoid  (under  the  natural  operation  (A)  •  (B)  :=  ( AB )).  Moreover,  the 
previous  operation  partitions  the  invertible  n  x  n  circulant  matrices  into 
equivalence  classes,  say  C*/~,  and  consequently,  (C*/~  ,  •)  becomes  a  group. 

Let  f  —  x\Xj2  *  *  *  'i' j,i  A  X‘) -V j 2 -}- 1  *  *  *  '  1  R  j‘2 — 1 '  *  *  j ,j — 1  b e  an 
MRS  function  of  degree  d,  with  the  SANF  x\ Xj2  ■  ■  ■  Xjd.  We  associate  to  / 
the  following  (unique)  circulant  matrix  equivalence  class 


1 

A/  =  (C(1,0,.. 


32  3 3  3d 

1 , 0, ....  0,  1 0,  1 , . . .  ,0)>, 


(1) 


where  the  1  bits  (indicated  above)  appear  in  positions  given  by  the  indices 
in  the  SANF  monomial  of  /. 

For  a  binary  (row)  vector  (ai,  <22, . . . ,  an)  of  dimension  n,  we  let 
A(ai,  a2, ... ,  an)  =  {*|  ai  =  1},  and  by  abuse  of  notation,  A(C(a))  =  A(a). 
Similarly,  for  a  single  monomial  term  x'q  X{2  ■  ■  ■  Xid  of  degree  d  in  n  vari¬ 
ables,  we  define  A (xqx^  •  •  •  Xid)  =  {ij\  j  =  1,2, ,  d}.  We  can  also  ex¬ 
tend  this  to  the  MRS  function  with  this  SANF,  /  =  xq Xi2  •  ■  ■  xid ,  as 
A(/)  =  A(xh  xi2  ■  •  ■  xid)i  which  is  not  unique,  but  we  prefer  (so  not  to 
complicate  the  notation)  to  consider  all  such  sets  equal  under  a  cyclic 
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rotation  permutation  of  the  indices.  That  is,  for  Af  as  in  (1),  then 
A{f)  =  {1,J2,  •  •  • ,  jd}  =  {2,  jd  +  1}  =  ••••  Therefore,  for  a 

fixed  n,  any  such  particular  set  A  of  indices  generates  a  unique  monomial 
xtl  xi2  ■  ■  ■  Xid  in  n  binary  variables,  a  unique  n-dimensional  bit  vector  a,  the 
corresponding  unique  matrix  equivalence  class  (7(a),  and  the  corresponding 
unique  MRS  function  /  =  x^Xi2  ■  ■  ■  Xid  (SANF)  such  that  Af  =  (7(a). 

We  now  consider  another  type  of  equivalence  between  circulant  matri¬ 
ces,  that  can  be  extended  to  the  equivalence  classes  we  have  defined.  For 
two  circulant  matrices  A,  B,  if  there  are  permutation  matrices  P,  Q  such 
that  PA  =  BQ,  then  A  and  B  are  called  P-Q  equivalent.  It  is  known 
in  that  case  that  AAT  and  BBT  are  similar  matrices  (in  fact,  there  ex¬ 
ists  a  permutation  matrix  which  conjugates  one  to  the  other)  [12].  More¬ 
over,  it  is  rather  straightforward  to  see  that  AAT  =  YlijeA(A)  ,  where 
A  —  C  (or , . . . ,  cin ) . 

The  notion  of  P-Q  equivalence  extends  naturally  from  circulant  matri¬ 
ces  to  equivalence  classes,  as  any  product  of  permutation  matrices  is  also 
a  permutation  matrix,  and  any  two  representative  matrices  A\,A2  of  an 
equivalence  class  (A)  are  related  by  a  rotation  of  the  row  order. 

The  following  two  results  are  essential  in  our  investigation. 

Theorem  2.1  (Canright-Chung-Stanica  [2]).  Two  MRS  Boolean  functions 
f,g  in  n  variables  are  S-equivalent  if  and  only  if  their  corresponding  circu¬ 
lant  matrix  equivalence  classes  Af  and  Ag  are  P-Q  equivalent. 

Theorem  2.2  (Wiedemann-Zieve  [12]).  Let  A,  B  be  two  nxn  0/1  -circulants 
of  weight  at  most  5  with  first  rows  support  indices  A(A),  respectively,  A (B), 
where  n  is  odd.  Then  the  following  are  equivalent: 

1.  There  exist  u,v  £  Zn  such  that  gcd(it,  n)  =  1  and  A(A)  =  uA(B)  +  v. 

2.  A,  B  are  P-Q  equivalent. 

3.  There  is  an  nxn  permutation  matrix  P  such  that  AA7  =  PBB 1  P~7 . 

4.  The  matrices  AAr ,  BBT  are  similar. 

It  is  known  [12]  that  if  the  weight  of  A,  B  is  2,  3,  then  Theorem  2.2  holds 
without  constraints.  If  the  weight  is  k  €  {4,5},  then  Theorem  2.2  holds 
for  dimensions  n  whose  prime  factors  are  greater  than  2 k(k  —  1)  (that  is,  if 
k  =  4,  the  prime  factors  should  be  greater  than  24). 

For  easy  reference,  based  upon  the  previous  theorems,  we  often  write 
interchangeably  /  ~  g  (for  two  MRS  /,  g),  or  A  (Af)  ~  A(T9)  (if  there  exist 
u,v  £  Zn  such  that  gcd(u,  n)  =  1  and  A(T)  =  uA(B)  +  v ). 

In  this  paper  we  will  find  the  number  of  equivalence  classes  (and  repre¬ 
sentatives  of  these  classes)  for  quartic  (degree  4)  MRS  (that  is,  their  SANF 
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is  /  =  x\ XiXjXk  with  A (/)  =  {1  ,i,j,k})  in  prime  and  prime  power  dimen¬ 
sions.  We  mention  that  the  result  for  prime  dimension  appears  in  a  recent 
paper  of  Cusick  and  Cheon  [6].  We  give  here  a  rather  short  proof  of  that 
result. 


3.  Counting  quartic  equivalence  classes  for  prime  dimension 

For  easy  displaying,  we  sometimes  write  f  to  mean  ab~ 1  in  some  obvious 
environment;  we  also  adopt  the  convention  throughout  that  working  in  some 
Zpe,  ptx~ 1  exists  if  pa\\x  =  pay,  0  <  a  <  t,  and  ptx~1  :=  pt~ay^1. 

We  start  with  a  descriptive  lemma  describing  some  representatives  of 
the  equivalence  classes.  Throughout  this  paper  we  use  the  “capital  mod” 
notation  a  Mod  n  to  mean  the  unique  integer  b  £  {1,2, ...  ,n}  such  that 
b  =  a  mod  n.  We  use  the  notation  ps\\m.  to  mean  the  highest  power  of  p 
that  divides  m. 


Lemma  3.1.  The  S -equivalence  class  of  any  quartic  MRS  h  with  A (h)  = 
{1  where  either  gcd(i  —  l,n)  =  gcd(j  —  l,n)  =  gcd(fc  —  l,n)  = 

gcd(i  —  j,  n)  =  gcd (k  —  j,  n )  =  gcd(A:  —  i,  n)  =  1,  contains  a  quartic  MRS  g 
with  A (g)  =  {l,2,m,  r}.  If  n  =  pe ,  l  >  2,  where  p  is  an  odd  prime  and 
gcd(i  —  l,n)  /  1,  gcd(j  —  l,n)  ^  1,  gcd (k  —  1  ,n)  1,  then  the  class  of  h 

will  not  contain  any  MRS  function  g  with  A(g)  =  {1,2 ,m,r},  rather  it  will 
contain  one  with  A(g)  =  {l,ps  +  1,  m,  r},  where  ps||  gcd(i  —  1,  j  —  1,  k  —  1), 
1  A  s  <  £  —  1,  ps\(m  —  1)  and  ps\(r  —  1).  Moreover,  every  such  class  with 
a  representative  of  support  {1  ,ps  +  1  ,Aps  +  1  ,Bps  +  1}  will  contain  one 
with  support  {l,ps  +  1,  aps  +  1,  bps  +  1}  with  gcd (ab,p)  =  1  where  p  divides 
either  both  of  (a  —  1) ,  (6  —  1)  or  none,  or  p\  gcd(a,  b  —  1)  or  p\  gcd(6,  a  —  1). 
Furthermore,  two  MRS  with  supports  {l,ps  +  l,aps  +  l,bps  +  l},  gcd (ab,p)  = 
1,  respectively,  {1  ,ps  +  1,  a'ps+l  +  1,  ( b'p+  1  )ps  +  1}  cannot  be  S -equivalent. 

Proof.  We  first  assume  that  at  least  one  of  gcd(f— 1,  n )  =  1,  gcd(j  — 1,  n)  =  1, 
gcd (k  —  1 ,  n)  =  1,  gcd(i  —  j,  n)  =  1,  gcd(/c  —  j,n)  =  1,  gcd(fc  —  i,n)  =  1 
holds.  By  Theorem  2.1  and  Theorem  2.2  it  will  be  sufficient  to  show  that 
for  every  such  MRS  li  with  A(/t)  =  {1  ,i,j,k},  there  exists  u,  v  such  that 
uA(h)+v  =  {1,  2,  m ,  r},  for  some  m ,  r.  Solving  the  corresponding  24  systems 
we  obtain  the  possibilities  for  (m,r,u,v)\ 
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i-  1  1 

7 -  +  1,  7 - b  1,  - 7,  7 - 

i  —  /c  i  —  A?  k  —  %  %  —  k 

i  +  j  —  2k  i  —  1  1  i 

- - - , - :  +2,  - - , - : 

i  —  k  k  —  i  i  —  k  k  —  i 

i+j-  2k  j  -  1  2  1  i 

j  —  fc  ’  fc  —  j  ’  j  —  fc  ’  fc  —  j 


+  2/ 

+  2}; 

+  2}; 

+  2l. 


We  easily  see  that  if  gcd(i  —  1 ,  n)  =  1,  taking  u  =  (i  —  l)-1,n  =  u  —  1, 
then  A:}  ~  {1,2 ,m,r}  via  u,v,  similarly,  for  gcd(j  —  1 ,  n)  =  1,  or 

gcd (k  —  1,  n)  =  1,  etc. 

Next  assume  that  for  s  >  1,  ps||  gcd(i  —  1,  j  —  1,  k  —  1)  (and  consequently, 
it  also  divides  j  —  i,k  —  i,k  —  j).  Without  loss  of  generality,  we  assume  that 
ps\\(i  —  1),  and  so  i  —  1  =  pst  for  some  t  ^  0  (mod  p)  (the  other  cases  are 
similar).  By  taking  u  =  t_1,  v  =  1  —  u,  m  =  1  +  (j  —  l)u,r  =  1  +  (k  —  1  )u 
(all  are  Mod  //),  then  we  see  that  {l,i,j,  k}  ~  {l,ps  +  l,m,r}.  Certainly 
ps\m  —  1  =  (j  —  l)t_1,ps|r  —  1  =  (k  —  l)t_1,  since  ps||  gcd(i  —  1,  j  —  1,  k  —  1). 
Moreover,  for  every  class,  which  we  have  now  shown  that  has  a  representative 
of  the  form  {1  ,ps  +  1,  aps  +  1,  bps  +  1},  if  p\a,  we  can  find  a  representative 
{1  ,ps  +  1  ,a'ps  +  1  ,b’ps  +  1}  with  p  /[ a'  (that  is  seen  by  taking  the  affine 
transformation  based  on  it  =  —  l,v  =  2  +  ps  and  so,  a'  =  1  —  a,  b'  =  1  —  b, 
gcd  (a',p)  =  1). 

The  next  claim  follows  easily  by  finding  appropriate  values  of  u,  v,  a,  b 
such  that  u{l,ps  +  1,  Aps  +  1,  Bps  +  1}  +  v  =  {l,ps  +  1,  aps  +  1,  bps  +  1} 
with  gcd (ab,p)  =  1  and  p\  gcd(a  —  1,  b  —  1),  or  p\ (a  —  1)  and  p\b,  or  p\(b  —  1) 
and  p\a : 

(i)  If  p  /\A(A  —  1  )B(B  —  1),  respectively,  we  can  take  u  =  (B  —  l)-1, 
v  =  1  -  (ps  +  1)(AB  -  I)-1,  a  =  -(B  -  I)"1,  b  =  (A  -  1  )(B  -  1)_1,  so 
p  J(ab(a  —  1  )(b—  1)  and  p  /(a  —  b)  (Case  B3  and  B5  later). 

(ii)  If  p\(A  —  l),p\(B  —  1),  then  p\(A  —  B ),  p  )(AB\  If  p\A,p\B,  then  we  can 
take  u  =  —  1,  v  =  2  +ps,  a  =  1  —  A,  b  =  1  —  B;  thus,  p\  ( a  —  l),p\(b  —  1),  and 
p\(a  —  b),  p  J(ab  (both  of  these  instances  are  equivalent  and  constitute  Case 


6 


B\  later). 

(in)  If  p\(B  —  1)  and  p  /A(A  —  1)  (similarly  for  p\(A  —  1)  and  p  )(B(B  —  1)), 
then  take  u  =  A ~1,v  =  1  —  A _1,  a  =  A~l  ,b  =  BA~l\  thus,  p\(a  —  b ), 
p  j(ab(a  —  1  )(b—  1)  (Case  B2  and  B±  later). 

(■ iv )  If  p\A,  p  )(B,  p  J((B  —  1)  (similarly,  for  p\B,  p  /A,  p  J((A  —  1)),  the  same 
transformation  as  in  (i)  gives  us  a  representative  with  p  f[ab(a  —  1  )(b  —  1) 
and  p\(a  —  b). 

( v )  If  p\A,p\  ( B  —  1)  (similarly  for  p\B,p\  (A  —  1)),  then  we  have  a  representa¬ 
tive  of  the  form  {l,ps  +  1,  aps+1  + 1,  (bp+  l)ps  +  1}  (a  =  A/p,b  =  (B  —  1  )/p) 
(Case  C  later). 

Thus,  a  representative  of  such  a  class  can  be  taken  of  the  form  {l,ps  + 
1,  aps  + 1,  bps  + 1}  with  gcd (ab,p)  =  1,  p\  gcd (b—  1,  a  —  1),  or  p  /(a  — 1)(6—  1)). 

Now,  given  two  classes  with  supports  {l,ps  +  1  ,aps  +  1  ,bps  +  1}  with 
gcd (ab,p)  =  1,  p\(a  —  1)  or  ( b  —  1),  respectively,  {l,ps  +  1  ,a'ps+l  +  1,  (b'p  + 
1  )ps  +  1},  and  listing  the  possibilities  of  ( a',b',u )  with  gcd (u,p)  =  1  that 
could  potentially  map  the  first  support  unto  the  second  via  the  correspond¬ 
ing  affine  transformation  with  some  shift  v,  we  see  that  in  every  instance, 
either  a,  or  b  must  be  divisible  by  p.  and  that  is  impossible.  Thus,  the  two 
classes  are  not  S'-equi valent.  □ 


We  denote  by  E(p)k  the  number  of  distinct  equivalence  classes  of  quartic 
MRS  in  p  variables,  for  p  =  k  (mod  12),  where  k  =  1,  5,  7, 11.  Although  this 
result  is  shown  in  [6],  we  give  a  much  shorter  alternate  proof  here. 


Theorem  3.2.  Let  p  >  11  be  a  prime, 
classes  of  quartic  MRS  in  p  variables  is 


E(p)  1  = 
E(p)7  = 


p2  —  2p  +  25 
24  ' 

p2  —  2p  +  13 
24  ’ 


Then  the  number  of  S -equivalence 


E(p)5 

E(p)u 


p2  —  2p  +  9 
24  ’ 

p2  —  2p  —  3 
'  24 


Proof.  Since  p  is  prime  by  Lemma  3.1  it  is  sufficient  to  find  the  number 
of  nonequivalent  MRS  with  support  {l,2,m,  r}.  For  that  purpose,  we  Lx 
3  <  j  <  k  <  p  and  look  at  possible  3  <  m  <  r  <  p  such  that  {1,  2,  j.  k}  ~ 
{1,2,  m,  r}.  Solving  the  corresponding  systems  we  obtain  the  following  12 
putative  values  of  {m,  r}  (unordered  pairs): 


{j,k};{3- j,3- k};{l  + (j  -  1)  L  1  +  (k  —  l)(j  —  1)  x}; 

{1  +  (fc  -  l)-1, 1  +  0  -  l)(fc  -  l)'1}  ;  {2  -  (j  -  l)-1, 2  -  (k  -  1  )(j  -  l)-1}  ; 

{2  -  (fc  —  1)  — 1 , 2  —  (j  -  1  )(k  -  I)"1}  ;  (1  -  (j  -  2)  — 1 . 1  +  (k  -  2 )(j  -  2)  1 } ;  (2) 

{1  -  (fc  -  2)  1 , 1  +  C j  -  2 )(k  -  2)  — 1 }  ;  (2  +  {j  -  2)  — 1 , 2  —  (fc  -  2 )(j  -  2)"1}; 
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{2  +  (k  —  2)“\  2  -  (j  -  2)(fc  -  2)  — 1 }  ;  {l  -  (j  -  1  )(fc  -  j)-\  1  —  0'  —  2)(fc  -  D”1}  i 
{!  +  (*-  i)(fc  -  i)_1, 1  +  (*  -  2)(fc  -  i)-1}  ; 


Case  1.  Let  k  =  j  +  1,  3  <  j  <  p  —  1.  Then  the  possible  values  Mod  p  for 
{m,  r}  are: 

{3,3  +  1};  (3  -3,2-  j}-,{  1  +  (j  -  I)"1,  2  +  (j  -  l)-1};  {1  +  j“\  2  -  j’1}; 

{2  -  (j  -  I)"1,  1  -  (j  -  I)"1};  {2  -  j~\  1  +J-1}; 

{1  -  (j  -  2) — 1 ,  2  +  (j  -  2)"1};  {1  -  (j  -  l)-\2  -  (j  -  I)'1}; 

{2  +  (j  -  2)-\  1  -  (j  -  2)"1};  {2  +  (j  -  1)~\  1  +  (j-  I)'1}; 

{2  ./.3  ./}:  {  1 

and  removing  the  obvious  duplications,  we  get 

{3,3  +  Ml  (3  -3,2-  j}-,{  1  +  (j  -  I)"1,  2  +  (j  -  l)-1};  {1  +  j~\  2  -  j'1}; 

{2  -  (j  -  I)"1, 1  -  (j  -  l)"1};  {1  -  C 3  ~  2)-\2  +  ( j  -  2)-1}. 

(Observe  that  four  pairs  of  consecutive  indices  are  contained  in  each  such 
class.) 

If  p  =  1,5  (mod  12),  by  Gauss’  reciprocity  law,  —1  is  a  quadratic  residue 
modulo  p,  and  so,  for  j  =  1  ±  (— l)1/2  Mod  p  (which  happens  when  j  = 
1  —  (j  —  l)-1,  for  example),  the  above  set  of  possible  values  for  {m,r} 
shrinks  to 

{j,  j  +  1};  {3  -  j,  2  —  j}-  {1  +  j- L,  2  -  j- L}. 

Certainly,  if  j i  =  l  +  (— l)1/2,  j‘2  =  1  —  ( — l)1/2,  then  the  class  of  {1,  2,  ji,ji  + 
1}  is  the  same  as  the  class  of  {l,2,ji,j2  +  1}.  (Observe  that  two  pairs  of 
consecutive  indices  are  contained  in  each  such  class.) 

The  contributions  to  E(p )(.)  in  all  these  cases  are 


e(p)  1,5  <-  1  +  - 

.  p —  3 

^(p)r,n  ^ - 4— 


P-1 

4 


(3) 


Case  2.  Assume  k  7^  j  +  1  (of  course,  by  symmetry,  nor  j  7^  k  +  1).  If 
p  =  1,7  (mod  12),  then  —3  is  a  quadratic  residue  modulo  p.  For  {j,  /c}  = 
{(9  ±  (— 3)1//2)6_1,  (3  ±  (— 3)1//2)2-1}  Mod  p  (these  values  are  obtained  by 
solving  the  system  3  —  j  =  1  +  (j  —  !)(&:  —  l)”1, 3  —  k  =  1  +  (&  —  l)-1,  for 


example),  the  set  (2)  shrinks  into  a  set  of  4  elements 


{.?,*};  {3-  j,3-fc}; 

{1  +  c?  -  l)-1, 1  +  (fc  -  1)(J  -  l)-1}  ; 

In  all  the  other  cases  when  k^j  +  1,  j^k  +  1,  either  they  all  belong  to 
the  above  class,  or  the  cardinality  of  the  set  (2)  is  12. 

The  contributions  to  both  E(p)\,  E(p)j  in  this  last  case  is 


1. 


(4) 


Putting  together  (3)  and  (4)  we  obtain 


PY  ,  _P-1  .  .  .  (V) -3-4-6CJ2  _p2-2p  +  25 
E(ph-  —  +  1  + - - - — 

p~2)-3-6  2=5  p2-2p  +  9 


E^  =  P-T1  +  ^^2 


24 


p  2)  -  4  -  6  ^  _  p2  -  2p  +  13 


E(ph  =  i^  +  1  +  ^— - __ 

FM  _p-3  .  (r?)-ep-^  _P*-2p-3 

£(p)i.-  — + - J2 - 55 - , 


24 


which  proves  our  theorem. 


□ 


4.  Counting  quartic  equivalence  classes  for  prime  power  dimen¬ 
sions 

We  will  now  count  the  S'-equivalence  classes  for  quartics  in  p(  variables 
(£  >  2),  where  p  >  29  is  a  prime  number.  We  start  with  a  few  lemmas. 

Lemma  4.1.  Let  f  be  a  quartic  MRS  inpe  (£  >  2)  dimension  whose  support 
includes  {1,2 ,j,k},  where  k  =  j  +  1  and  p\(j  —  1),  or  p\(j  —  2),  or  gcd(j  — 
1  ,p)  =  gcd(j  —  2 ,p)  =  1.  Then  its  equivalence  class  contains  an  MRS  whose 
support  is  {1,2 ,j',k'}  where  j'  +  k!  =  3  (modpf).  Furthermore,  the  class 
containing  {1,2 ,j,k}  with  j  +  k  =  3  (mod  pf)  will  not  contain  a  class  of 
support  {1,2 ,j',k'}  with  gcd (A/  —  1  ,p)  =  gcd (k1  —  2 ,p)  =  gcd (j'  —  1  ,p)  = 
gcd (/  -2,p)  =  gcd (k1  -  f,p)  =  1,  | k'  -  f\  /  1. 

Proof.  If  p\  (j  —  1)  we  take  the  affine  transformation  of  Theorem  2.2  based 
upon  u=(j-  2)_1,  v  =  1  -  2 (j  -  2)_1  and  {/,  A:'}  =  {1  -  (j  -  2)^x,  2  +  (j  - 
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2)-1}  (thus,  j'  +  k'  =  3  (mod  jr);  we  also  note  that  p\(j'  —  1)  and  p\(k'  —  2), 
or  p\(f  ~  2)  and  p\(k'  —  1)).  For  the  second  claim,  that  is,  if  p\(j  —  2), 
we  take  the  affine  transformation  based  upon  u  =  =  2  +  j”1  and 

{f,  k'}  =  {2  -  j-1, 1  +  j-1}  (thus,  f  +  k'  =  3  (mod  /)).  If  p  /[ (j  -  1), 
p  /(j  —  2),  then  either  of  the  two  above  transformations  will  do  the  trick. 

The  last  claim  follows  easily  by  solving  the  corresponding  system  and 
analyzing  each  of  the  twelve  solutions.  Since  we  have  done  several  such  and 
will  do  more  later,  we  leave  this  as  a  warm-up  exercise  for  the  reader.  □ 

Lemma  4.2.  Let  f  be  a  quartic  MRS  in  pe  (l  >  2)  dimension,  whose 
support  includes  {1,2 ,j,k},  where  p  divides  one  of  ( j  —  1  ),(j  —  2),  and 
gcd (k  —  1  ,p)  =  gcd (k  —  2 ,p)  =  1;  or,  p  divides  one  of  ( k  —  1  ),(k  —  2), 
and  gcd(j  —  1  ,p)  =  gcd(j  —  2,p)  =  1;  or  p  divides  both  (j  —  1),  ( k  —  1); 
or  p  divides  both  (j  —  2),  (k  —  2).  Then  its  class  contains  an  MRS  with 
support  {1,2,/,  k'}  where  p\(k'  —  f)  and  gcd  (A/  —  1  ,p)  =  gcd  (k'  —  2  ,p)  = 
gcd  (j'  -  1  ,p)  =  gcd  (f  -  2,p)  =  1. 

Proof.  If  p |  (j  —  2)  and  gcd (k  —  1,  p)  =  gcd (k  —  2 ,p)  =  1,  we  can  use  the  affine 
transformation  based  on  u  =  (k  —  =  1  —  (k  —  l)-1,  so  {f,k'}  =  {k{k  — 

1) -1,  l  +  (i-l)(^’-l)_1};  ifpIO-2)  and  gcd(j-l,p)  =  gcd(j'-2,p)  =  1,  we 
can  use  the  affine  transformation  based  onit=  (j  —  l)-1,  v  =  1  —  (j  —  l)-1, 
so  { j',k'j  =  {j(j  -  I)"1, 1  +  (k  -  1  )(j  -  l)-1};  if  p\(j  -  1)  and  gcd  (A;  - 
1  ,p)  =  gcd(k  —  2 ,p)  =  1,  we  can  use  the  affine  transformation  based  on 
u  =  (k  -  2)-\u  =  1  -  2 (k  -  2)_1,  so  {f,  k'}  =  {1  -  (k  -  2)“\  1  +  {j  - 

2) {k  —  2)”1};  if  p\(k  —  1)  and  gcd(j  —  1  ,p)  =  gcd(j  —  2 ,p)  =  1,  we  can  use 
the  affine  transformation  based  on  u  =  ( j  —  2)~1,v  =  1  —  2(j  —  2)_1,  so 
{/,  k'}  =  {1  —  ( j  —  2)~l ,  1  +  {k  —  2 )(j  —  2)-1}.  The  lemma  is  shown.  □ 

Remark  4.3.  If  p  J(B  =  A  +  1  (similarly,  for  p  /A  =  B  +  1)  then  {1  ,ps  + 
1,  Aps  +  1,  Bps  +  1}  ~  {1  ,ps  +  1,  aps  +  1,  bps  +  1}  with  a  +  b  =  1  (mod  pe~s ) 
(use  the  transformation  based  onu  =  —  (^4  + 1)-1,  v  =  1  —  (A  + 1)_1,  a  =  l  — 
(^4-t-l)-1,  b  =  (1— a)  Mod  pe~s ).  Also,  p  /(a  — 1)6,  and  if  p  f{A—l)A(A+l), 
then  p  fab(a  —  1)(6  —  1). 

Theorem  4.4.  Let  p  >  29  be  a  prime.  Then  the  number  of  S -equivalence 
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classes  of  quartic  MRS  in  pl  variables  (t  >  2)  is 

p2e+ 2  +  p2l+1  +  pn  -  3 pe+2  -  6 pe+1  -  3 pe  +  p2{27t  +  2)  +  5p  -  27£  +  2 

24  (p2  -  1) 

p2£+2  +  p2^+i  +  p2r  _  3^+2  _  6pe+ 1  _  3/  +  p2(n^  +  2)  +  5p  -  11£  +  2 

24  ( p 2  -  1) 

p2«+2  +  p2^+i  +  p2r  _  3^+2  _  6j/+i  _  zpe  +  p2(15e  +  2)  +  5p  -  15£  +  2 

24  ( p 2  -  1) 

p2«+2  +  p2^+i  +  p2^  _  3pf+2  _  6^+1  _  3pe  -  p2(l  -  2)  +  5p  +  ^  +  2 

24  (p2  -  1)  ' 

Proof.  By  Lemma  3.1,  to  count  the  number  of  equivalence  classes  of  quar- 
tics  in  pr  dimension,  it  will  be  sufficient  to  count  the  classes  with  sup¬ 
port  {1,2 ,j,k}  with  3  <  j,  bk  <  p,  and  also  count  classes  with  support 
{1  ,ps  +  1,  aps  +  1,  bps  +  1},  gcd(a6,p)  =  1,  1  <  s  <  i  —  1,  or  p\  gcd(a,  6—1) 
(or  p\  gcd(6,  a  —  1)),  2  <  a,  6  <  pe~s  —  1. 

Case  A.  First,  given  a  class  with  support  {1,2, j,  k},  the  only  possible 
values  for  {m,  r}  with  {1,  2,  j.  k}  ~  {1,  2,  m,  r}  are  listed  in  (2).  We  consider 
several  subcases. 

Subcase  (^4i).  k  +  j  =  3  (mod  p£)  (of  course,  j  /  k).  We  get  the  fol¬ 
lowing  possibilities  (removing  the  obvious  duplications)  for  {m,  r}  where 
{1,2,  j,  k}  ~  {1,2 ,m,r}: 

{j,  3  -  -  l)-1, 1  +  0  -  l)-1};  {-U  -  2)"1, 1  -  0  -  2)-1}; 

{2  -  (j  -  l)-\  3  -  ( j  - 1)"1};  {2  +  (j  -  2)— 1 , 3  +  (. j  -  2)"1}; 

{(3  -  (2 j  -  3)-1)  2“\  (3  +  (2 j  -  3)-1)  2"1}. 

(Observe  that  2 j  —  3  is  invertible,  unless  j  =  1  +  2_1  =  PM—  (mod  p);  also, 
j  ^  (mod  //),  since  then  j  =  3  —  j  Mod  //,  which  is  impossible.) 

It  is  straightforward  to  check  that  there  are  exactly  two  pairs  {m,  r}  in 
the  above  list  satisfying  m.  +  r  =  3  (mod/),  namely  the  first  and  the 

l,o 

last  ones,  except  when  j  =  2_+_  (mod/,  or  when  p  =  1,5  (mod  12) 
(since,  by  Gauss’  reciprocity  law,  —1  is  a  quadratic  residue  modulo  p) 
and  so,  for  j  =  (3  ±  (— 1)1/,2)2_1  Mod  p  the  two  pairs  overlap  (precisely, 
j  =  (3  —  (2 j  —  3)_1)  2_1  and  3  —  j  =  (3  +  (2 j  —  3)”1)  2”1).  It  is  obvious 

£  Q 

that  there  are  P-fA  ordered  pairs  (j,  k)  ( j ,  k  >  3)  with  k  +  j  =  3  (mod  p 
out  of  which  p  2-1  of  them  satisfy  j  =  £-±5  (mod  p). 


E(p%  = 
E{p%  = 
E(pe)r  = 
E(pe)  n  = 


11 


Thus,  the  contribution  to  E(jr)^  in  this  case  is 

t~  1  _  1  1  /  mf-  _  Q  mi-t  _  X 


W)l,5 


-®(P^)7,11 


1  + 


p 


-1  ,  1//-3  p' 


p^+p*"1  -2 


P 


l—l 


—  1  1  /  p  —  3  p4  1  —  1 


pf+p£-1-4 


(5) 


Subcase  (A2).  k+j  ^  3  (mod  j/),  p|  (/c— j),  and  gcd(/c—  l,p)  =  gcd(k— 2,p)  = 
gcd(j-l,p)  =  gcd(j-2,p)  =  1,  orp|  gcd(j  -  1,  k  -  1),  orp|  gcd(j  -  2,  k  -  2). 
We  get  the  following  possibilities  for  {m,  r}  where  {1,2,  j,  k}  (1,2  ,m,r}: 

for  gcd(fc  -  l,p)  =  gcd(fc  -  2 ,p)  =  gcd(j  -  l,p)  =  gcd(j  -  2 ,p)  =  1: 

O'. {3  -3, 3  -  fc};  {1  +  O'  - 1)^1, 1  +  (k  - 1)0  -  I)-1} ; 

{1  +  (fc  -  l)-1, 1  +  0  -  m  -  l)-1}  ;  {2  -  U  -  l)-1, 2  -  (fc  -  1)0  -  l)"1}  ; 

{2  -  (fc  -  I)”1,  2  —  0  —  l)(fc  -  I)”1}  ;  {1  -  0  -  2)“\  1  +  (fc  -  2)(j  -  2)’1}; 

{1  -  (k  -  2)"1, 1  +  0  -  2 ){k  -  2)-1}  ;  {2  +  0  -  2)"1, 2  -  (k  -  2)0  -  2)'1}; 

{2  +  (A;  —  2)_1,  2  —  0  —  2)(fc  —  2)-1}  , 

and 


for  p|0  —  1)  an<i  p|(fo  —  1): 

0.  kh  (3  -3,3-  fc};  (1  -  0  -  2)-1, 1  +  (k  -  2)0  -  2)-1}; 

{1  -  (fc  -  2)-1, 1  +  0  -  2)(fc  -  2)-1}  ;  {2  +  0  -  2)"1, 2  -  (k  -  2 )(j  -  2)"1}; 

{2  +  (fc  —  2)_1,  2  —  0  —  2)(fc  —  2)-1}  , 

for  p|0  —  2)  and  p\(k  —  2): 

O.  kh  (3  -3, 3  -  fc};  {1  +  0  -  1)“\  1  +  (k-  1)0  -  I)-1} ; 

{1  +  (fc  - 1)-\  1  +  0  -  m  -  I)-1} ;  {2  -  0  - 1)-1, 2  -  (fc  -  1)0  - 1)-1} ; 

{2  —  (fc  —  1)1, 2  —  0  —  i)(fc  —  i)-1}  • 

We  observe  that  the  last  two  possibilities  are  not  different  since,  if  p\(k  —  1) 
and  p\(j  —  1),  then  the  class  of  {1,2,  j,  k}  will  contain  {1,2,  j',  k1}  with 
p\(k'  —  2),p\(j'  —  2)  (for  example,  {f,  k'}  :=  {3  —  j,  3  —  k}),  and  viceversa. 
Moreover,  there  are  exactly  two  pairs  in  each  class  satisfying  the  required 
conditions;  for  example,  if  p\(j  —  l),p\(k  —  1),  then  in  the  class  of  {1,  2,  j,  k} 
one  can  find  only  two  others  {1,2 ,j',k}  with  p\(j'  —  l),p\(k'  —  1),  namely, 
if,k'}  =  {2  +  (j-2)~\2-(k-2)(j-2)-1}  and  {f,k'}  =  {2+{k-2)~\2- 
(j-2)(k-2)~'}. 

Summarizing,  every  class  in  the  first  category  (gcd(fc  —  1  ,p)  =  gcd(fc  — 
2,  p)  =  gcd  (j  —  1  ,p)  =  gcd  (j  —  2  ,p)  =  1)  contains  only  two  pairs  (recall  that 
j  +  k  ^  3  (mod  pr),  as  in  that  case  we  would  have  only  one  pair),  satisfying 
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the  imposed  conditions,  and  every  class  in  the  second  category  (p\ (j  —  1)  and 
p\(k  —  1),  or  p\(j  —  2)  and  p\(k  —  2))  contains  exactly  three  pairs  satisfying 
the  imposed  conditions. 

There  are  — - ^4^ - —  ordered  pairs  (j,  k)  satisfying  p\ (j  —  1) ,  p\  (k  —  1) 

and  the  same  number  satisfying  p\ (j  —  2) ,  p\  ( k  —  2) .  There  are  - - j- - - 

ordered  pairs  (j,  k)  with  p\ ( k  —  j ) ,  from  which  we  take  away  the  ones  satisfy¬ 
ing  p\(j  —  1)  and  p\(k  —  1),  or  p\  ( j  —  2)  and  p\(k  —  2),  or  j  +  k  =  3  (mod  pe). 
Using  an  inclusion- exclusion  easy  argument,  we  obtain 

(j/-1  -  l)  ( pe  -  4)  o  ( pe _1  -  1  )(pe~1  -  2)  pE~l  -  2 

2  2  2 
(]/_1  —  1  )(pe  —  2pe~1  —  1) 

~~  2 

number  of  pairs  (j,  k)  with  p\(k  —  j)  /  0  and  j  +  k  ^  3  (mod  p£)  such  that 
gcd(fc  -  l,p)  =  gcd(fc  -2,p)  =  gcd (j  -  1  ,p)  =  gcd (j  -  2,p)  =  1. 

Thus,  the  contribution  to  E(p £)(.)  in  this  case  is 


1  (p*-1  -  1  )(pe  -  2pi~1  -  1)  1  (/-1  -  l)^"1  -  2) 

2  2  ^  3  2 

(p^-1  -  l)(3p£  -  4/"1  -  7) 

_  12  ' 


Subcase  (VI3).  j  +  k  ^  3  (mod  pe),  \j  —  k\  /  1,  and  p\(k  —  1)  and  p\(j  —  2), 
or,  p\  ( k  —  2)  and  p\  ( j  —  1)  (thus,  p  J((k  —  j)).  The  possible  values  for  {m,  r} 
with  {1,2 ,j,k}  ~  {1,2,  m,  r}  are  (all  different  Mod  p£): 

for  p\(k  —  2)  and  p\(j  —  1): 

{j,  fc};  {3  -  j,  3  -  k}- {1  +  (k  -  I)"1, 1  +  0  -  l)(fc  -  I)"1}  ; 

{2  -  (k  -  I)"1, 2  -  0  -  i)(fc  -  I)"1}  ;  {1  -  0  -  2)-\  1  +  (k  -  2)0  -  2)-1}; 

{2  +  0  -  2)  1 , 2  —  (k  -  2)0  -  2)"1};  {l  -  (j  -  1  )(fc  -  j)-\l  -  0  -  2)(fe  -  j)"1}  ; 

{1  +  (k  -  l)(fc  -  j)~\  1  +  (*  -  2)(fc  -  j)-1}  . 

for  p|(fc  —  1)  and  p|(j  —  2): 

0,  fc};  {3  -j,  3  -  fc};  {1  +  0  -  1)_1, 1  +  (*  - 1)0'  - 1)-1} ; 

{2  -  0  -  1)_1,2  -  (k  -  1)0  -  l)-1}  ;  {1  -  (k  -  2)-1, 1  +  0  -  2 )(k  -  2)-1}  ; 

{2  +  (k  -  2)-\  2  -  0  -  2)(fc  -  2)"1}  ;  {l  -  (j  -  1  )(k  -  j)~\  1  -  ( j  -  2)(k  -  j)-1}  ; 

{1  +  (k  -  1  )(k  -  j)~\l  +  (fc  -  2)(fc  -  jy1}  . 

We  observe  that  the  two  possibilities  are  not  different  since,  if  p\[k  —  1) 
and  p\(j  —  2),  then  the  class  of  {1,2 ,j,k}  will  contain  {1,2,  f,  k'}  with 
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p\(k'  -  2),p\(f  -  1)  (for  example,  k'  :=  1  +  (j  -  1)  \  f  :=  1  +  {k  -  l)(j  - 
l)”1),  and  viceversa.  Moreover,  every  pair  {f,  k'}  in  the  class  of  {j,  k}  with 
p\(k  —  1)  and  p\(j  —  2)  satisfies  the  same  condition  (for  some  choice  of  f,  k'): 
for  example,  for  \k'  :=  1  +  {k  —  1  :=  1  +  ( j  —  1  )(k  —  1)— 1 } ,  we  have 

p\(f  ~  2),p|(fc'  -  1). 

The  number  of  pairs  {j,k}  ( j,k  >  3,  \j  —  k\  /  1)  with  p\(k  —  2)  and 
p\(j  —  1)  and  k  +  j  ^  3  (mod  //)  is  exactly  ( p —  1  ){p(’~l  —  3).  Thus,  the 
contribution  to  E(pe)^  in  this  case  is 


E(jf 


1.5,7,11 


(P 


-1 


Vi/-1  ~3) 
8 


(7) 


Subcase  (A4).  k  +  j  ^  3  (mod  pe),  gcd (k  —  1  ,p)  =  gcd (k  —  2 ,p)  =  gcd(j  — 
1  ,p)  =  gcd  (j  -  2,p)  =  gcd  (k  -  j,p )  =  1,  as  well  as  |  k  -  j\  p  1  (j,  k  >  3). 
The  possible  values  for  {m,r}  Mod  p(  with  {1,2 ,j,k}  ~  {1,2,  m,r}  are 
as  in  (2).  As  before,  if  p  =  1,7  (mod  12),  then  —3  is  a  quadratic  residue 
modulo  p;  for  {j,k}  =  {(9  ±  (— 3)1/,2)6_1,  (3  ±  (— 3)1/2)2“1}  Mod  pe  (or, 
{j,  k}  =  {(5  ±  ( — 3)  1/2)2~1 ,  (5  ( — 3)1/2)2— 1},  etc.),  the  set  (2)  of  values  for 

{m,  r}  shrinks  into  the  same  set  of  four  elements 

{j,  k};{3  -  j',3  -  k};{l  +  ( j  -  l)-1, 1  +  (k  —  1  )(j  -  l)-1}  ; 

{2  -  (i  -  l)"1, 2  -  (fc  -  l)(i  -  I)"1}  . 

Since  we  often  use  the  inclusion-exclusion  principle  in  the  paper  (in  the 
“underground”)  we  will  display  explicitly  such  a  use  of  the  it  to  count  the 
number  of  ordered  pairs  with  the  above  conditions  (we  always  assume  that 
3  <  j  <  k  <  pl). 

Let  X\  =  {(j,  k)  |  \k  —  j |  ±  1,  k  +  j  =  3  (mod  /)},  A2  =  {(j,  k)  \  \k  - 
j  I  +  l,p\k  ~  l, or  p\k  -  2},  X3  =  {(j,k)  \  \k  -  j  \  /  l,p\j  -  l, or  p\j  -  2}, 
X4  =  {(j,k)\  \k  —  j\  /  l,p\k  —  j}.  We  see  that 

Ui|  =  'krl,  \x2\  =  Us|  =  (P‘  -  - 1),  \xt\  =  kl - 4>y~1  - 11 

(the  last  count  is  obtained,  by  observing  that  for  every  j  >  3,  we  let  k  =  ap+ 

j,  where  1  <  a  <  pe~l  —  j/p,  and  so,  \X4\  =  YH=i  1p(pt~1  —  i)  —  2(pf~1  —  1), 
which  simplifies  to  the  above  expression).  Observe  that  the  universal  set 
S  =  {(j.  k)  |  k  >  j  +  1  >  4}  has  \S\  =  (p  {"2)  —  (//  —  3)  (to  account  for 
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k  =  j  +  1).  Further, 


|*i  ni2|=  vl~x  -  1;  |*i  n  x3\  =  p €~1  -  l; 

|Xi  n  x4\  =  pi  12~1-,  \X2  n  X3\  =  2(/~1  -  l)^"1  -  2); 

\x2  n  x4\  =  (/-1  -  l)^-1  -  2);  \x3  n  XA\  =  (p*-1  -  l)!/-1  -  2); 
|Xi  ni2ni3|  =  /_1  -  1;  |*i  ni2n x4|  =  o; 

\Xi  m3n xA\  =  O;  \x2  ni3n xA\  =  ft/-1  -  i)(P£~l  -  2); 

\Xi  nx2nx3nx4\  =  o. 


Using  the  complementary  form  of  the  inclusion-exclusion  principle  for  N  =  4 
sets  in  the  universal  set  S  we  introduced  earlier,  we  obtain 


n* 


N 


=  \s\-J2  l-W  +  I]  I Xi  nXj\ - +  (-i)n|Xi  n . . .  n  xN\ 

1  <i<j<N 

(pe  —  4)(p£_1  —  1) 


i= 1 

p"  -  2\  3 (pe  -  3) 


-  2 (pe  -  5 )(/  1  -  1)  - 


+ 


3  (P(  ~  1) 


+  3(/-1-l)(p"-1-2) 


p21  —  5 p2t  1  +  6 p2t  2  —  3 pe  +  9 pe  1  _  pl  1  (pe  —  2 pe  1  —  3)  (p  —  3) 


Thus,  the  contribution  to  E(pe)^  in  this  case  is 


E(p£) 
E(pl ) 


1,7 


5,11 


1  + 


pl  1  [pf  —  2  p 


l- 1 


24 


jr  1  (jr  —  2jr  1  — 


24 


~3)(p~3) 

3)(p-3) 


1 

3’ 


(8) 


Case  B.  Next,  fix  1  <  s  <  £  —  1  and  consider  a  class  of  support  {1  ,ps  + 
1,  aps+l,  bps+l},  gcd(ab,p)  =  1  (with  p  dividing  either  both  of  (a— 1),  (6—1), 
or  none).  We  consider  separately  the  two  subcases  (these  are  disjoint,  as 
one  can  see  from  the  proof  of  Lemma  3.1).  (To  simplify  our  counts  a  bit,  we 
assume  that  all  the  expressions  are  reduced  modulo  jr~s,  because  aps  +  1 
is  the  same  index  as  (a  +  p^~s)ps  +  1  Mod  // . )  By  abuse  of  notation,  if 
{1  ,ps  +  1,  aps  +  1,  bps  +  1}  ~  {l,ps  +  1,  mps  +  1,  rps  +  1},  then  we  say  that 
{a,  b},  {to,  r}  are  equivalent  and  often  write  {a,  6}  ~  {to,  r}. 

Subcase  (B i).  p  J(ab,  p\(a—l)  and  p\(b—  1).  The  potential  values  of  {to,  r}  ~ 
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{a,  6}  Mod  pe  s  are 

{a,  b};  {1  —  a,  1  —  b};  {a-1,  6a-1};  {ft-1, aft-1}; 
{1  —  a-1, 1  —  fta-1};  {1  —  ft-1, 1  —  aft-1}. 


(We  make  here  an  observation,  which  we  will  use  later:  all  six  values  { m ,  r} 
above  satisfy  the  property  that  p\(m  —  ?’).)  Removing  the  above  pairs  that 
do  not  obviously  satisfy  the  conditions p\(a—l),p\(b—l)  we  get  the  potential 
values  (all  distinct  Mod  pr~8): 


{a,  ft};  {a  1,6a  x};{6  1,aft  x}. 

For  a  fixed  s,  there  are  — - -1^p - ordered  pairs  (a,  ft)  with 

p |  gcd(a  —1,6—1),  and  so,  summing  them  over  1  <  s  <  i  —  1,  we  obtain 
that  the  contribution  to  E(pe)^  in  this  subcase  is 

p2e~2  -  3/  -  3 p*-1  +  2(i  -  1  )p2  -  2(1  -  2)  +  3 p 
E(p  )i,5,7,n  < - - 71 - 

6(P  1}  (9) 

_  (/-1  -  1 )(//  :"l  -  3p  -  2)  l  -  1 

6  (p2  -  1)  +  3 

Using  the  above  observation,  there  are  — - 3f>  ^  +^  —  1  pairs  (a,  ft) 

withp|(a  — 6)  andp  /(a— l)(ft— 1)  (since  for  all,  p  divides  both  a,  ft)  contained 
in  all  of  these  classes,  which  we  have  to  disregard  later. 

Subcase  (B2).  p  /aft(a  —  l)(ft  —  1),  p\  (a  —  6),  a  +  ft  ^  1  (mod  pe~s).  The  ten 
potential  values  Mod  pi  of  {a,  ft}  are  (all  distinct) 

{a,  ft};  {1  —  a,  1  —  6};  {a-1, 6a-1};  {fo  1 ,  afo—  1 } ; 

{1  —  a-1, 1  —  6a-1};  {1  —  6-1, 1  —  a6-1}; 

{-(a  -  l)-1,  (6  -  l)(a  -  l)-1};  {-(6  -  l)-1,  (a  -  1)(6  -  l)-1}; 

{1  +  (a  -  l)-1, 1  -  (6  -  l)(a  -  l)-1};  {1  +  (6  -  l)-1, 1  -  (a  -  1)(6  -  l)-1}. 


We  observe  that  in  the  above  list  there  are  only  two  pairs  that  fully  sat¬ 
isfy  the  imposed  conditions.  An  elementary  inclusion-exclusion  argument 

shows  that  there  are  — -  1  2 - - - -  for  which  p\(a  —  ft),  a  +  ft  ^  1 

(mod  pe~s),  p  /aft(a  —  1)(6  —  1).  Summing  over  1  <  s  <  £  —  1,  we  obtain 
a  total  of  — — — o(^'-\)+P — ~P+2  +  |  such  pairs.  The  contribution  to 
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E(pe)^  in  this  case  is 


E{pt)  1,5,7,11 


p2e  1  —  2 p21  2  —  pe+1  +pe  1  —  p  +  2  i 

4(p2  -  1)  +  4 

(p2e~2  -  l){p  -  2) 

4(p2  -  1)  4  ' 


(10) 


Subcase  (£>3).  p  )(ab(a  —  l)(b—l)(a  —  b),  a  +  b  ^  1  (mod  p£~s).  If  \a  —  b\  =  1, 
say  b  =  a  +  1,  there  exists  a  representative  based  on  a',  b'  in  the  same  class, 
satisfying  the  conditions  p  )(a'b'(a'  —  l){b'  —  l)(a7  —  b')  along  with  \a'  —  b'\  >  1 
(use  transformation  u  =  b~l,v  =  1  —  ft-1,  a'  =  b~1,b/  =  ab -1).  Thus,  we 
will  assume  that  \a  —  b\  >  1. 

The  potential  values  of  {m,  r}  Mod  pr,  where  {l,ps  +  l,  aps  + 1,  bps+ 1}  ~ 
{l,ps  +  1,  mps  +  1  ,rps  +  1}  (observe  that  2  <  a,b  <  pe~s)  are: 

{a,  6};  {1  —  a,  1  —  6};  {a-1,  6a-1};  {6-1,  a&-1};  {1  —  a_1,l  —  6a-1}; 

{1  -  6_1, 1  -  ab-1}- {-(a  -  I)'1,  (6  -  l)(a  -  l)-1};  {-(6  -  l)-1,  (a  -  1)(6  -  l)-1}; 

{1  +  (a  —  l)1, 1  —  (6—  l)(a  —  1) — 1 } ;  { 1  +  (6—  l)-1, 1  —  (a  —  1)  (6  —  l)-1};  (11) 

{a(a  —  ft)-1,  (a  —  l)(a  —  6)-1};  {6(6  —  a)-1,  (6  —  1)(6  —  a)-1}. 

If  p  =  1,7  (mod  12),  then  there  exist  two  (nontrivial)  cubic  roots  a,  b  of 
1  (this  is  equivalent  to  pe  \  a2  +  a  +  1 ,  p£  \  b2  +  b  +  1 ;  we  also  note  that  the 
two  cubic  roots  satisfy  a  +  b  +  1  =  0  (mod  p^~s))-  For  these  two  values 
(they  occur  when  the  first  and  fourth  pairs  in  (11)  are  the  same),  of  course, 
reduced  modulo  p^~s,  the  set  (11)  shrinks  into  a  four  element  set 

{a,  6};  {1  -  a,  1  -  6};  {-(a  -  l)-1,  ( b  -  l)(a  -  l)-1};  {-( b  -  l)-1,  (a  -  1)(6  -  l)-1}. 


In  all  the  other  cases,  the  set  (11)  has  12  distinct  elements. 

Yet  another  inclusion-exclusion  argument  reveals  that  there  are  ns  := 

- - — — - — 3)(p— 3)  orc(ereci  pairs  (a,  b)  with  p  J(ab(a  —  1  )(b  —  l)(a  —  b ) 

and  |a  —  b\  ^  1,  a  +  b  ^  1  (mod  p£-s).  Summing  over  1  <  s  <  l  —  1,  we 

obtain  a  total  of  ^  - 2P2(p2-i) — ~ — +2'P+5'>  such  pairs.  For  p  =  1,7 

(mod  12),  we  sum  ns  —  4  for  every  s,  to  account  for  the  pairs  in  the  class  of 
the  mentioned  cubic  roots. 

The  contribution  to  E(pe)^  in  this  case  is 


E(pe) i,7  <—  (<-!)  + 


(p  -  3)  [p21-1  -  2 p2l~2  -  3 pl  -  3 p1-1  +  2p  +  5)  4(£  -  1) 


24  (p2  -  1) 


12 


E(p%,  i 


(p  -  3)  (p21-1  -  2 p2l~2  -  3 pl  -  3 p1-1  +  2p  +  5) 
24  (p2  -  1) 


(12) 
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Subcase  (£>4).  p  f\ab(a  —  1)(6  —  1),  p\(a  —  b),  a  +  b  =  1  (mod  pe  s ).  The 
set  (11)  contracts  into  the  5  element  set 

{a,  1  -  a};{— 1  +  aT1,  a”1};  {— (a  —  l)-1,— 1  +  (a  -  l)'1}; 

{1  -  a-1, 2  -  a'1};  {1  +  (a  -  1)_1, 2  +  (a  -  l)-1}. 


We  note  that  in  the  above  list,  there  is  only  one  pair  satisfying  the  imposed 

conditions,  and  there  are  - — g -  such  ordered  pairs,  which  summed  for 

1  <  s  <  £  —  1  gives  the  contribution  to  E(pe)  in  this  case 


E(pe)i,  5,7,11 


P(P*  2  -  1) 
2(p-l) 


1  —  2 
2 


(13) 


Subcase  {B§).  p  J(ab(a—  1)(6  —  l)(a  —  b),  a  +  b  =  1  (mod  j/  s).  The  set  (11) 
shrinks  to  six  elements,  namely 

{a,  1  -  a};  {a-1,  -1  +  a-1};  {-(a  -  1)“\  -1  -  (a  -  l)-1};  {1  -  a-1, 2  -  a-1}; 

{1  +  (a  -  I)-1, 2  +  (a  -  I)"1};  {a(2a  -  l)”1, 1  -  a(2a  -  l)”1}. 

We  note  that  only  two  pairs  in  the  above  set  satisfy  the  imposed  conditions 
(the  first  and  the  last  pairs).  If  p  =  1,5  (mod  12),  then  —1  is  a  quadratic 
residue  modulo  p£,  and  for  {a,  6}  =  {1  ±  (— 1)1/,2)2-1}  (where  {a,  1  —  a}  = 
{1  —  a(2a  —  l)_1,a(2a  —  l)'1})  the  above  set  shrinks  into  the  same  four 
element  set  (regardless  of  the  chosen  sign) 

{a,  6};  {1  -  a,  1  -  6};  {a-1,  6a'1};  {1  +  (6  -  1)_1, 1  -  (a  -  1)(6  -  l)'1}. 


In  all  the  other  cases  the  set  (14)  contains  distinct  elements.  There  are 

ms  :=  - - 2^- — -  pairs  satisfying  the  conditions  of  this  case,  which  summed 

over  1  <  s  <  Z  —  1,  gives  us  a  total  of  —  2^-1)  ^  such  pairs.  For  p  =  1,5 
(mod  12)  we  sum  ms  —  2  (to  account  for  the  pair  {1  ±  (— 1)1//2)2-1}),  for 
every  s. 

Thus,  the  contribution  to  E(p( in  this  case  is 


E(p\  5 
E(p% 


11 


(p  -  l)(p~  3)  l-l 
4  (p  -  1)  2 

(/-1  -  l)(p  —  3) 

4(p  -  1) 


(15) 


Case  C.  Last,  fix  1  <  s  <  £—2  (of  course,  this  case  needs  £  >  3)  and  consider 
a  class  of  support  {1  ,ps  +  1,  aps+1  + 1,  {bp+  1  )ps  +  1}  ( a '  :=  ap ,  6'  :=  6p+  1). 
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The  8  potential  values  of  (a,  6)  Mod  p{  s  1  are  (we  consider  ordered  pairs 
here,  since  a  appears  in  aps+1,  while  b  is  in  (bp  +  1  )ps,  so  their  roles  are  not 
interchangeable) 

(a,  6);  (-b,  -a);  (a{bp  +  1)_1,  -b(bp  +  1)_1);  (b(bp  +  1)_1,  -a(bp  +  1)_1); 

(a{ap  -  1)_1,  -b(ap  -  1)  x);  (b{ap  -  1)_1,  -a(ap  -  1)_1);  (16) 

(—a(bp  —  ap+  1)_1,  —b(bp  —  ap  +  1)_1);  (■ b(bp  —  ap  +  l)-x,  a(bp  —  ap  +  1)_1). 


We  note  that  (assuming  order)  the  set  (16)  contains  either  4  or  8  disjoint 
elements.  We  have  4  elements  when  the  class  has  a  representative  (a,  b)  with 
a  +  b  =  0  (mod  p f^s^1),  or  a  =  b.  The  number  of  such  pairs  is  exactly 
2(p*-s_1  _  2)  (the  two  conditions  do  not  overlap,  since  if  a  =  b  and  a  + 
6  =  0  (mod  p^_s_1),  then  it  follows  that  p^^1  \  a,  which  is  impossible). 
The  remaining  number  of  pairs  (of  class  cardinality  8)  is  exactly  (p^_s_1  — 
1 )(pe~s~1  -  2)  -  (p*-®-1  -  1)  =  p2e-2s~ 2  -  4 p1-*-1  +  3. 

Thus,  the  contribution  to  E(p ^)(.)  of  class  C  is 


p _ 2 

E{pe)l,5,7,ll  < - ^  f 

s=l  ' 


-l  V 


21— 2s— 2 


—  4 p‘ 


1- s-1 


y?  p2^-2S-2  _  1 
5  =  1 


J2i-2 


- 1 


1 


(17) 


8  (p2  -  1)  8 


Putting  together  equations  (5),  (6),  (7),  (8),  (9),  (10),  (12),  (13),  (15)  and 
(17),  we  obtain 


„  t,  ,  pe+pe~1-2  ,  (p*-1-!)^- V_1-7)  ,  (p1-1-  Vtf-'-S) 

E(P  )i  < - 4 - + - la - + - § - 

p/-1(pf-2p<-1-3)(p-3)  1  (p1-1  -Vtf-1  -Zp-2)  l-l 

+  24  3  6  (p2  -  1)  3 

(p2e~2  -  l)(p-2)  (p  —  3)  (p2f_1  —  2p2e~2  —  3pe  —  3pe~1  +  2p  +  5) 

4  4(p2  -  1)  4  1  24(p2  -  1) 

.  m^-l)  1-2  (/-‘-ll^S)  M  p2»-l  l-l 

1  ’  2(p  —  1)  2  4(p  —  1)  2  +8(p2-l)  8 

_  p2l+2  +  p2i+1  +  p2t  -  3 pe+2  -  6pl+1  -  3pe  +  p2  (27 1  +  2)  +  5p  -  277  +  2 
“  24  (p2  -  1)  ’ 

p‘+p‘-1-2  (p<-1-l)(3pf-4p<-1-7)  (p^1-l)(p^-1-3) 

4  12  8 


E(pl 


|  pe~1(pt  —  2pe~1  —  3)(p—  3)  |  (p"~1-l)(p^1-3p-2)  |  l-l 


24 


6  (p2  —  1) 


(p2e~2  —  1)  (p  —  2)  pe~1-£  (p  —  3)  (p2f_1  —  2p2l~2  —  3pe  —  3pe_1  +  2p  +  5) 

J  /  O  -I  \  4 


4(p2  —  1) 


24(p2  -  1) 
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P(pe~2~i)  1-2  (p1-1  -  i)(p  ~  3)  e-i  p2e~2  —  1  e- i 

+  2(p-l)  2  +  4(p-l)  +  2  +  8(p2-l)  8 

_  p2f+2  +  p2l+1  +  p21  -  3 pe+2  -  6 pt+1  -  3pe  +p2{lli  +  2)  +  5p  -  1U  +  2 

2K//'  -  1)  ’ 

,  ,  pe  +  pei  —  4  (p^-lHS/Z-V-1-?)  (p^-lXp^-S) 

)7  4 - 4 - + - 12 - + - 8 - 

1  pe~1(pe  —  2pe~1  —  3)(p  —  3)  1  (/-1-l)(p<-1-3p^2)  l-l 

24  3  6  (p2  -  1)  3 

(p2e~2  -  l)(p-2)  p^"1-^  (p  —  3)  (p2^-1  —  2p2l~2  —  “ip1  —  3pe_1  +  2p  +  5) 

4“  77  o  7  4~ 


4(p2  —  1) 


24(p2  -  1) 


I  n  |  p(/~2-l)  *~2  (p^-lXp-S)  p2^-2-i  i-  1 

1  j  2(p-l)  2  4(p  —  1)  8(p2-l)  8 

_  p2t+2  +  p2t+1  +  p21  -  3 pt+2  -  6 pe+1  -  3 pe  +  p2(  15£  +  2)  +  5p  -  1U  +  2 

24  (p2  -  1) 

^  ,  P*  +  /'1  -  4  ,  (p*"1  -  l)(3pf  -  4 p1-1  -  7)  ,  (p^1  -  i; )(p^1  -  3) 

£(P  )n  4 - 4 - + - 12 - + - 8 - 

|  pe~1(pe  —  2p1^1  —  3)(p—  3)  |  (pf~1-l)(p<-1^3p-2)  |  €-1 


24 


6  (p2  —  1) 


{p2l~2  —  l)(p  —  2)  p4"1-^  (p  —  3)  (p2t~1  —  2p2l~2  —  3pe  —  3pe_1  +  2p  +  5) 

+  4(p2  -  1)  4  +  24  (p2  -  1) 

p(/-2-  1)  £-2  (pt_1  —  l)(p  —  3)  p2f-2-l  l-l 

+  2(p-l)  2  4(p  —  1)  8(p2-l)  8 

p2t+2  +  p2t+i  +  p2i  _  3^+2  _  6pm  _  Sp(  _  p^£  _  2)  +  5p  + 1  +  2 

24  (p2  -  1) 


The  theorem  is  shown. 


□ 


We  end  with  Table  1  displaying  the  number  of  equivalence  classes  for 
dimension  //,  where  3  <  p  <  31,  and  powers  1  <  t  <  5.  As  an  example, 
the  actual  computation  (using  a  program  supplied  by  T.W.  Cusick  and  his 
students)  for  n  =  73  took  10.5  hours  on  an  Windows  7  Pro  i7  with  16GB  of 
memory,  however,  for  n  =  5  ,  the  program  did  not  finish  the  computation 
after  6  days.  Interestingly,  our  formulas  agree  with  all  these  computations  for 
smaller  powers  (in  spite  of  the  fact  that  we  had  to  impose  the  bound  p  >  29 
in  our  main  theorem)  and  so,  we  used  our  formulas  to  quickly  compute  the 
number  of  classes  for  higher  powers  (we  put  an  asterisk  next  to  that  count). 

Acknowledgement.  The  author  would  like  to  thank  T.W.  Cusick  for  some 
useful  suggestions  and  for  his  help  with  a  Mathematica  code.  Also,  thanks 
are  due  to  the  referee  for  providing  some  useful  suggestions. 
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Tabic  1:  The  number  of  equivalence  classes  for  pe,  3  <  p  <  31,  1  <  £  <  5 


n 

P 

1 

2 

3 

4 

5 

3 

0 

4 

44 

426 

3940 

5 

1 

30 

819 

20908* 

524997* 

7 

2 

112 

5766 

284840* 

13973842* 

11 

4 

658 

81612* 

9897016* 

1197780520* 

13 

7 

1274 

218757* 

37019332* 

6256916099* 

17 

11 

3670 

1071393* 

309820076* 

89541196087* 

19 

14 

5698 

2073642* 

748909856* 

270362647330* 

23 

20 

12140* 

6458556* 

3417415664* 

1807832195324* 

29 

33 

30446* 

25695699* 

21612733932* 

18176386152005* 

31 

38 

39676* 

38246514* 

36758592152* 

35325121572190* 

[1]  A.  Biryukov,  C.  De  Canniere,  A.  Braeken,  B.  Preneel,  A  Toolbox  for 
Cryptanalysis:  Linear  and  Affine  Equivalence  Algorithms,  Advances  in 
Cryptology  -  Eurocrypt  2003,  LNCS  2656  (E.  Biham,  ed.),  Springer- 
Verlag,  2003,  33-50. 

[2]  D.  Canright,  J.H.  Chung,  P.  Stanica,  Circulant  Matrices  and  Affine 
Equivalence  of  Monomial  Rotation  Symmetric  Boolean  Functions,  to 
appear. 

[3]  T.W.  Cusick,  Affine  equivalence  of  cubic  homogeneous  rotation  sym¬ 
metric  functions,  Inform.  Sci.  181:22  (2011),  5067-5083. 

[4]  T.W.  Cusick,  A.  Brown,  Affine  equivalence  for  rotation  symmetric 
Boolean  functions  with  pk  variables,  Finite  Fields  Applic.  18:3  (2012), 
547-562. 

[5]  T.W.  Cusick  and  Y.  Cheon.  Affine  equivalence  for  rotation  symmetric 
Boolean  functions  with  2k  variables,  Designs,  Codes  and  Cryptography 
63  (2012),  273-294. 

[6]  T.W.  Cusick,  Y.  Cheon  Affine  equivalence  of  quartic  homogeneous 
rotation  symmetric  Boolean  functions,  Inform.  Sci.  259  (2014),  192- 
211. 


21 


[7]  T.W.  Cusick,  P.  Stanica,  Fast  Evaluation,  Weights  and  Nonlinearity  of 
Rotation-Symmetric  Functions,  Discrete  Mathematics  258  (2002),  289- 
301. 

[8]  A.  Maximov,  M.  Hell,  S.  Maitra,  Plateaued  Rotation  Symmetric 
Boolean  Functions  on  Odd  Number  of  Variables,  International  Work¬ 
shop  on  Boolean  Functions:  Cryptography  and  Applications  (BFCA 
2005),  University  of  Rouen,  France  (2005);  available  at  eprint.iacr.org, 
no.  2004/144,  2004. 

[9]  A.  Maximov,  Classes  of  Plateaued  Rotation  Symmetric  Boolean  func¬ 
tions  under  Transformation  of  Walsh  Spectra,  Workshop  on  Coding  and 
Cryptology  2005  (O.  Ytrehus,  ed.),  LNCS  3969  (2006),  325-334. 

[10]  J.  Pieprzyk,  C.X.  Qu,  Fast  Hashing  and  Rotation-Symmetric  Functions, 
J.  Universal  Computer  Science  5  (1999),  20-31. 

[11]  P.  Stanica,  S.  Maitra,  Rotation  Symmetric  Boolean  Functions  -  Count 
and  Cryptographic  Properties,  Discrete  Appl.  Math.  156  (2008),  1567- 
1580. 

[12]  D.  Wiedemann,  M.E.  Zieve,  Equivalence  of  sparse  circulants:  the  bi¬ 
partite  Adam  problem,  http://arxiv.org/abs/0706.1567. 


22 


